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We investigate evolutional paths of an extended quintessence with a non-minimally coupled phan- 
tom scalar field ip to the Ricci curvature. The dynamical system methods are used to investigate 
typical regimes of dynamics at the late time. We demonstrate that there are two generic types of 
evolutional scenarios which approach the attractor (a focus or a node type critical point) in the phase 
space: the quasi-oscillatory and monotonic trajectories approach to the attractor which represents 
the FRW model with the cosmological constant. We demonstrate that dynamical system admits 
invariant two-dimensional submanifold and discussion that which cosmological scenario is realized 
depends on behavior of the system on the phase plane We formulate simple conditions on 

the value of coupling constant ^ for which trajectories tend to the focus in the phase plane and hence 
damping oscillations around the mysterious value w = — 1. We describe this condition in terms of 
slow-roll parameters calculated at the critical point. We discover that the generic trajectories in 
the focus-attractor scenario come from the unstable node. It is also investigated the exact form of 
the parametrization of the equation of state parameter w{z) (directly determined from dynamics) 
which assumes a different form for both scenarios. 

PACS numbers: 98.80.Bp, 98.80.Cq, 11. 15. Ex 



I. INTRODUCTION 



The non-minimal coupling of the type ^£_Rip'^ in the Lagrangian, where R is the Ricci scalar and ^ - the parameter 
of non-minimal coupling was independently introduced by Chernikov and Tagirov y and Callan et al. It was 
pointed out that we should consider ^ 7^ and different arguments were given. The nonzero ^ arises from quantum 
corrections [3] and it is required by the renormalization [5] . Faraoni argued that the non- minimal coupling of the 
inflaton field is unavoidable and in his opinion the paradigm of standard inflation should be generalized by including 
the coupling term in the Lagrangian. If the parameter ^ assumes a value different from 1/6 (conformal coupling case) 
then the equation of motion for the scalar field (i.e., the Klein-Gordon equation) is conformally invariant if Viip) = 
or V{'ip) (x tp'^ While there are strong arguments that seems to favor the choice of ^ = 1/6, under the assumption 
of conformal invariance of the Klein-Gordon equation and that ip does not violate the equivalence principle [g, 0| 
(contrary to Lightman et al. d, p. 85]). Moreover the value of £ = 1/6 cannot be treated as a correct coupling for 
various scalar field particles. Adopting a different point of view Q it is desirable to study different constraints on the 
coupling constant. The strength of the non-minimal coupling is constrained to be ^ > —7.0 * 10"'^ for negative ^ in 
the case of a quadratic potential function [T^j . 

Recent astronomical observations (like SNIa, the CMB and other [ll|, [H) El) 13 ) indicate that current universe 
is almost flat and undergoing accelerated phase of expansion. While the cosmological constant term introduced to 
the equations of general relativity offers the possibility of explanation of acceleration, its value is extraordinary small 
compared with an inferred vacuum energy. To solve this problem it is natural to consider a source of this acceleration 
in the form of a scalar field (quintessence idea [l^ In the simplest case minimally coupled to gravity scalar 

fields with a potential V{4') are considered. In this case the dynamical effect of a scalar field is equivalent to a perfect 
fluid with the energy density and the pressure given by = ^-0^ -I- V^(V') and = ^ip"^ — V{ip), respectively. In 
the canonical quintessence scenario the kinetic energy is fc = "0^/2 and the true cosmological constant is zero and the 
potential (tracking potential) is slow rolling down [17| . 

This quintessence field generically has good-or-t)ad attractor properties [H, [l^, i.e., it is sensitive with respect 
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to choice of the initial conditions. The tracking quintessence was invoked to use this attractor property to explain 
the smallness of the cosmological constant without fine tuning of the potential or initial conditions. The idea of 
quintessence was extended to the case of non- minimally coupled scalar field [lO, HH, ■ The main motivation of this 
extension was to explain the observational data which favor the negative equation of state parameter w = p/ p < — 1 
[23| . If we consider the canonical scalar field in the background of Einstein gravity then we cannot explain the range 
of value w < — 1. If we assume the fiat universe then w < —1 implies superacceleration as a opposed to acceleration 
d/a = H + > 0, where a is the scale factor and H = (ln(a))' is the Hubble parameter. There are two natural 
possibilities of extension of the theory: 1) to allow the scalar field -0 to couple non-minimally to the Ricci scalar [2^ . 
and 2) including of the phantom scalar fields, i.e., the scalar field is formally allowed to have a negative kinetic energy 
and ip"^ is replaced by (— [H, [23|. This type of modification has a justification from string, M theory and in 
supergravity. Note that such possibilities can explain the breaking w < —1 barrier. Both phantom and non-minimally 
coupled scalar field models can be regarded as a scalar tensor theories of gravity. Such models are called extended 
quintessence [l^, 113, IH| ■ Recently Fuzfa and Alimi suggested a promising unified description of both dark matter 
and dark energy. It relies on a violation of the weak equivalence principle on cosmological scales by dark matter. It 
is interesting that amplitude of this violation depends on the relative concentration of baryons [29j . 

The main aim of this paper is to investigate the long time evolution of extended quintessence models. We give the 
complete qualitative classification of solutions according to critical point approach scenarios (the monotonic approach 
to a critical point or the quasi-oscillatory regime of approaching and crossing of the w — —1 barrier). We use the 
dynamical system methods in investigation of all evolutional paths for all admissible initial conditions. These methods 
were previously used (see for example (sol. [3l|) in the context of the FRW model with non-minimally coupled scalar 
fields and superacceleration. We are interested in typical evolutional scenarios of the phantom non-minimally coupled 
scalar field at a late time near the w = —1 state. We demonstrate that there are two different scenarios to achieve 
this state. The new scenario of the route to the cosmological constant is through the damping oscillations of the 
scalar field asymptotically going to the constant value. In this scenario there is an infinite number of crossing the 
phantom divide line w = — 1. In our previous paper (3^ we study how deSitter state can be achieved in cosmology 
with conformally coupled to gravity phantom scalar field. The dynamics of the phantom scalar field conformally 
coupled to gravity with quadratic potential function depends on the value of squared mass of the field. The dynamics 
is regular for > 0, in opposite case with < the dynamics becomes chaotic and fractals structure appears 
in the phase space [33| . In the present paper we generalize result of [33| by consideration ^ as a free parameter 
which should be constrained by observational data. In investigation of dynamics with the help of dynamical system 
methods this parameter plays the role of bifurcation parameter and we found some critical value of this parameter for 
which topological structure of phase space changes. We show that damping oscillations around w — — I value type of 
evolution can be described by a critical point of a focus type on a 2-dimensional invariant submanifold (jp,ip') where 
' = d/dln(a). We will demonstrate that the late time evolutionary scenario crucially depends on the value of the 
non-minimal coupling constant. 



II. FRW MODEL WITH THE NON-MINIMALLY COUPLED PHANTOM SCALAR FIELD 

We assume the flat model with the FRW geometry, i.e., the line element has the form 

ds^ = -dt^ + a'^{t)[dr^ + r^{d0^ + sin^ M^^)], (1) 

where < (f < 2tt, Q < 9 < tt and < r < oo are comoving coordinates, t stands for the cosmological time. It is also 
assumed that a source of gravity is the phantom scalar field ip with an arbitrary coupling constant ^. The dynamics 
is governed by the action 

^^IJ d"a:V^(m2i? + (5'^^pV'. + C^?^' " 2C/(^))) (2) 

where = (SttG)^^; for simplicity and without lost of generality we assume 4ttG/3 = 1 and U{ip) is a scalar field 
potential. After dropping the full derivatives with respect to time, rescaling the time variable to the conformal time 
drj — dt/a we obtain the energy conservation condition 

£• = ifl^ + ^a^TA^ + S^a^V'^ + G^dipaijj - a^U{ip) - p,n,Qa = Pr^ (3) 

where pr,o and Pm,o are the constants corresponding to the radiation and the matter in the model, respectively, and 
dot denotes differentiation with respect to conformal time. We can also express this condition as a 

^H'^ ^ P^ + Pr + Pm (4) 
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where H is the Hubble function, is the energy density of a phantom scalar field, p,- oc a ■* is the energy density of 



radiation and p„ 



is the density of matter. 



The Euler-Lagrange equation reduces to 



a^-ip + G^aatp — ~2atpa + a'^U (1/1), 



(5) 



It is easy to check that for the conformal coupling ^ — 1/6, quadratic potential function, the rescaled scalar field 
variable ip = (f>/a and Pm,o = this system reduces to well known one. 
It would be useful to change 

a 

a a"^ 

where a dot denotes differentiation with respect to the conformal time and a prime with respect to a natural logarithm 
of the scale factor. Equations of motion (5) are in the form 



where 



U{^l)) + Pr+ Pn 



a4 1 + (1 - 6^)ip'^ + 6^(V'' + Tp)^ ' 
After eliminating the scale factor and it derivatives we obtain 

(V/' + [1 + 6^(1 - eov^'] + ^/'"(l - 60 (V' + 6^^-)+ 
(1 - eov-" + 6C(V'' + V')' 



(6) 



(7) 



1 1 
2~ 



(4(7(V^) + + 6^V) - C/'(^)(l + 6^V(V'' + ^)) 



= 0, 



(8) 



where a prime denotes the differentiation with respect to a natural logarithm of the scale factor. 

The dynamical systems methods offer a possibility of analyzing all evolutional paths for all admissible initial 
conditions in a phase space. They are represented by the trajectories in the phase space. The structure of the phase 
space is organized by the trajectories and limit sets. The critical points of the system are such points in the phase 
space for which right-hand sides of the dynamical system vanishes at these points. If the system is non-linear we can 
linearize it at these critical points and the famous Hartman-Grobman theorem guarantees that linearized system is 
a good approximation of the original dynamical system in the vicinity of this point. The type of a critical point is 
also characterized by a linearization matrix. It can be determined from the characteristic equation det |^ — Al| = 0, 
where A is the linearization matrix. The critical point is the global attractor if all real parts of the eigenvalues are 
negative, i.e., ReA^ < V^. 

Introducing the new variable y — ip' we can represent equation ([5]) as the autonomous dynamical system 



y' 



y 

ll + (l-60y^ + 6e(y + V)2 
2 



l + 6eV'2(l-60 



60 



U{ljj) + Pr + Pn 



(9) 



-4pr 

-3p„ 



where we have pm = Pm,oo, ^ and pr — Pr,oa, ^. 

The equation of state can be directly determined from the dynamics 

1 



WcS - 



1 -6C 



(1 + 2^rp\l - 6O) + e|J (pm + - 6^Vf/'(^))] }■ (10) 
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where 



1 + (1 - eov^'^ + + ipy 



is the counterpart of the Friedmann first integral. 

On the invariant submanifold we change the "time" variable 



da dln(a) 

where e = 1 when U{ip){l + 6^'i/'^(l ~ 6^)) > and e — —1 when U{i/j){1 + 6^V^(1 ~ 6^)) < 0. The dynamical system 
is in the form 

eip ^yU mi + 6^ij^ (1-60), (11) 
ey = - yC/W(l + G^i^^l - 6^) - (1 - GOv^V + 

- ^(1 + (1 - 60y' + 6^(2/ + ^)')(4C/(^)(y + 6^^) - U'^l + G^VX?/ + ^)))- 

In the generic case critical points of a dynamical system are located: 
1) yo = and (1 + 6^iji){24^^oUii;o) - C/'(V'o)(l + 6CV'o)) = «; 
^-^ ^0 — 6{(6^-i) for C > 1/6 or ^ < and yo is a solution of the cubic equation 

where for a given form of the potential function, A(^), -B(^), C(^), -D(C) are functions of the coupling constant ^ only. 

For our further investigations the critical point (1) is especially interesting. At this point the trace and the 
determinant of the linearization matrix are 

trA - -i7(^o)(l + 6eVo(l - 60) - ^(1 + 6eVg)(4C/(^o) - 6^%U'{^/jo)) (12) 



det A = --(1 + 6e^2)(i2eV'ot/'(V^o) + 24^C/(Vo) - U"{^„){1 + G^^j^oWi^JoKl + 6CVo'(l " 6^) (13) 

The type of behavior around this point crucially depends on the sign of A = (trA)^ — 4 det A. We are interested in 
the case ^ > 0. Note that at this critical point 



A 



^ _ U{^Po)U"{^Po) ^ 



Ui^Jo) 1 + 6e^2 
we obtain that 

(Ui^o))' (^(1 + 6CV^o'(l - 60) (9(1 + 6CVo(l - 60) - nm^i^l + C " r24^Vo))) (14) 
The general condition for A < is 

III. COSMOLOGY WITH NON-MINIMAL COUPLING AS A DYNAMICAL SYSTEM 

In this section we concentrate on investigation of the dynamics of the models with the quadratic potential because 
for other most popular forms of the potential functions there are attractors only for the negative coupling parameter. 
This case seems to be less interesting physically [l3]- For the quadratic potential U{ip) = ^m'^ip'^, = 1, the 
dynamical system on the invariant submanifold reduces to 

£V' = ^#'(l + 6^^2(^_g^))^ 

ey = -i#'(i + (1 - GOi^m^ + y) + y')) - ^^(i + (i - 60y' + Q^y + ^)')(2y^(i - 30 + 6^^^ _ (^g^ 
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where a dot means differentiation witfi respect to ^ = jV"^!! + 6^7/1^(1 — 6^))gT^[^- 
Tfie linearization matrix of system (16) at the critical point yo = and t/jq = ^ is 



A = 







(17) 



The eigenvalues of the linearization matrix are: 



Ai,; 



4C 



(1 - 30 ± ^^/A 



where A = (1 — 3^(3 — 25^. The type of the critical point depends on the value of A. 

• for < ^ < ^: 

A > and both eigenvalues are real and negative, the critical point is a stable node (FiglT]); 

• for^<e<i: 

A < we have two complex eigenvalues and the critical point is of a focus type (Figl2]); 

• for ^ > i: 

A > and Ai > 0, A2 < and the critical point is of a saddle type. 



• for C = ^ or^ = i: 

A = and critical point is degenerated (Fig[3]). 

The eigenvectors are 



Vl,2 



3 I \/3 / 3-25C 
2 ^ 2 V 1-35 



(18) 



They are helpful in the construction of the exact solution of the linearized system. For now, we restrict ourselves to 
the case of < ^ < 3/25, for other cases the construction of the linearized solution is similar. 



x{a) = x(0) exp <^a 
( 1 







1-3^ 

9 1-3? 



3 I V3 / 3-25^ _3 _ \/3 / 3-25; 
2 ^ 2 Y 1-3C 2 2 V 1-3? 



expAiCT 
expA2cr 



1 + \/I / 3-25g \/3 / 3-25g 

2 ^ 2 V 1-35 3 V 1-35 

1 _ \/3 / 3-25g _ \/3 / 3-25g 

2 2 Y 1-35 3 Y 1-3C 




(19) 



where ct = ^^'^ '^^^ xq = tpi — TpQ, yo — ip'^ — ip^ are initial conditions. The exact solutions are 



tp' — -ipQ = aia' 



/ 1 - 3C / 1 1 

3^ ^ 3^° 



\/3. 



' 1-3C /I 1 
3^ + 3^" 



2^0 



Vs. 



1-3C /I 



3 - 25e + 3^° 



Q;2a 



2^0 



1-3^ /I 



3-25e V2'^°"'3^° 



(20) 
(21) 



where xq, yo are initial conditions, ai — ^5 (3 — V3\/t^^) ^"^^ '^2 = ^5(3 + V3a/%:^-)- Notice that these 



solutions are also vaHd for ^ > 1/3 (a saddle type critical point). 

For the degenerated case ^ — 3/25, equations (20) and (21) simplify to 



2 " 



-3/2 



(22) 
(23) 
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The same procedure can be applied to construction the hnearized solution for A < 0. Instead we can use these 
solutions if we notice that for 3/25 < ^ < 1/3 



25^ 


- 3 




3^ 


25e 


- 3 




3^ 



— la, 



then we have 



xo cos (a In a) + 3a ^ sin (a In a) ( \xq -\- ij/o 



Ho cos (a In o) — ax^ sin (a In o) — -a ^ sin [a In a) ( -a;o + -yo 



^ - Vo = o"^/^ 



where xq, yo are initial conditions, and a depends on the coupling constant ^ only. 

The substitution of ijj and gives us a general formula for u>cff around the critical point of a given type. 



(24) 
(25) 



On the invariant submanifold and for quadratic potential function U{'4i) — \'^'^ formula for equation of state 



parameter reduces to 



-|i(l + - 60) - (1 - 20(1 + (1 - 60^'^) - 4^(1 - 30(V' + V')'} 



l + 6eV'(l-60 

In the case when the critical point is a stable node (0 < ^ < 3/25) we have 

- (1- 30 + /2(0a)a-3 



(1 - 30 + 601 - 60^-0(^0"' + Ba-«0a-3/2 + 3^(1 _ 60(Aa«i + Ba-'^^Ya'^ ' 



(26) 



(27) 



where Vo = ^, a; = ^^|^. 



25; 



3« 



, 1 AT / 1 - 30 1 1 \ 

A=-xo + V3W^-^(-xo + -yoj, 



xq and 2/0 are the initial conditions for and i)' , respectively, and functions /i and are given by 

/i(Oa) = 2$Vo((3(l - 40 - 4ai(l - 30) Aa"' + (3(1 - 40 +4ai(l - 30)Sa-"'), 
3 



/2(0 a) = ( - -(3 - 40 + 1501 - 20) (Aa'^' + Sa""') + 



+ a, (3(1 -40 -801 -30) (a 



2^2ai _ ^2^-20, 



(28) 
(29) 

(30) 
(31) 



This case can be slightly simplified if we notice that there are two independent directions leading to the critical point 
of the stable node type 



\.B = \xq 



V3 



/ 1 


-3? 




-255 


/ 1- 


-3« 



{\xq + lyo) = 0, 



Then 



(1 - 30 + V6?(l - 20xoa"' + (301 - 20 + SQl - 30ai + (1 - 40a,')xga^ 
(1 - 30 + \/6C(l - 60a;oa«* + 301 - 60a;ga2«i 



(32) 
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Using general formula 

Px = Pxfid exp I - 3 J ^, da | (33) 
we can calculate dark energy density in the case of monotonic approach to deSitter state 

Px = px,oA{Bo + Bia"' + Bsa^"-)'^ exp | - 7 arctan (^/l - 6^(1 + ^M^oa"')) } (34) 

where 

A={Bn + Bi+ B^y^ exp {7 arctan {^1~6^{1 + y^xo))}, 

Bo = 1 - 3^, Bi = V6e(l - 60^0, B2 = 3^(1 - ^0^1 and 

12^2 + 8^(1 -3e)a, + (l-4C)af 



/3 



2^(1 - 60a« 
-12e^ + 8e(l-3e)a. + (l-4e)Qf 



ai = (3 - x/3^i^ j and as - (s + \/3^t^ 

In the case of damped oscillatory approach to the critical point (3/25 < ^ < 1/3) and for critical point V'o — V^^i 
"00 = we have general equation of state parameter in the form 



o«c^ ~(1 - 3g) + gi(e, a)a-^'^ + ffzl^, a)a-^ 

(1 - 30 + 6^(1 - H)i^oh{i, a)a-y^ + 3^(1 - G^h^i^ a)a-^ 



(35) 



where the functions /i, gi and (72 are 



3 11 

h{£_,a) = xo cos (q!osc In a) H sin (aosc lna)(-a;o + 77 yo), (36) 

Qfosc 2 3 ' 

51 (0 a) = 2^-00(^(1 - 6^/1(^0) - 4(1 - 3C)((a;o + yo) cos (aosc In a) - aosca:^o sin (^osc Ina)- 

- ■;7-^sin(aosclna)(ixo + ^yo))), (37) 
2q!osc 2 3 V 

52(e,a) =e(l-60/^'(^,a)- 

/ 9 1 1 \2 

- (1 - 2^(1 - 60 2/oCOs(Q!osclna) - aoscaJo sin (aosc In a) - sin (aosc lna)(-a;o + 7:2/o) - 

V Zttosc 2 3 / 

/ 3 1 1 \2 

- 4^(1 - 30( (a;o + yo) cos (aosc In a) - aosc a;o sin (aosc In a) - 7; sin (aosc In a) (-xq + ttJ/o) , (38) 

V 2aosc 2 3 / 



where aosc = "^y T^3^ and xq, j/o and -00 have their usual meaning. 

In the case of conformal coupling and assuming that yo = —^xq, equation of state parameter (j35p simplifies to 



-I + ^a^^/^f cos (^ In a) + V7 sin (^ In a)) xo - ^a"^ f 4 - 3 cos (V71n a) + \/7 sin (v^ln a)) Xo- (39) 
3V2 2/6V / 



This expression can be easy integrated and we receive dark energy density in the case of conformally coupled phantom 
scalar field 

Px;5=i =Px,oexp - ^(^^"^ ~ 29xo)a;o) exp (^^a^^/^ (5 cos (^ In a) + V7 sin(^ \iia))xo^ 

exp (^a-^( - 32 + 3cos (y^lna) - 9V7sin {V7\na))xl) , (40) 
where xq and px.o ^re two constants fitted from observational data. 
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It is interesting to consider the special case of minimal coupling ^ = 0. The dynamical system is in the form 



dr ^U'i^P)' 



dy 
L dr 



(41) 



where -p- = rrn -r- 

d r U yip) d (T 



'^WJ^'dTna.- Above dynamical system can be easily integrated 



(42) 



y = iP' ^ tan (-SV^ + A t + Co), (43) 

where 

1 r d lnC/(V') , , , 
Ar.-y^ — ^dln«', (44) 

is the "period of time" between initial and final state and Co is the constant of integration. The equation of state 
parameter assumes very simple form for the minimal coupling 

?i;x,«=o = -l-2V". (45) 

For comparison, on the Fig.|4]we present phase portraits for minimally coupled phantom scalar field with compactifi- 
cation at infinity for various potential functions. It is interesting that for every discussed form of the potential there 
are always trajectories with wx,(,=o = — 1 at a certain moment of the evolution. 

In current astronomical analysis of the nature of dark energy, the most popular are linear (in the scale factor or 
the redshift) parametrizations of equation of state parameter. Let us compare equations ()27p and l|35p with current 
estimations of coefhcients of this parametrizations. Expanding ([TT]) and ([55]) in to the Taylor series with respect to 
the scale factor or the redshift at the present epoch we receive polynomial approximation of exact form of the equation 
of state parameter. As an example we use two linear parametrizations and values of wq and wi estimated by Dicus 
and Repko [S^l- In the Taylor series we have three free parameters ccq, yo and ^, but form observations we have only 
two numbers wq and wi, so one free parameter should be removed. Most convenient seems to be assumption that 
?/o = — f a^o- Taking in to account only two first terms in the Taylor series we can solve system on the ^ and xq. On the 
Fig. [5] we present dependence of xq on the coupling parameter ^ for the values of wq and wi form [34| . Intersections 
of this curves are solutions on xq and ^. It is interesting that the obtained values of the coupling parameter ^ with 
the assumption yo — — §2^0 are close to the conformal coupling — \- 



IV. ANALYSIS AT INFINITY 



The full description of a dynamical system includes also behavior of its trajectories at infinity. It is performed in 
tools of the Poincare sphere construction. In this approach we project trajectories from the center of the unit sphere 
onto the plane {x,y) tangent to either at the north or south pole (for details see [35l . p. 265]). Due to this 
construction infinitely distant points of the plane are mapped into the sphere equator, phase trajectories are mapped 
into corresponding curves on the sphere but the character of critical points is conserved and new critical points can 
appear at the equator. Hence the orthogonal projection of any hemisphere onto a tangent plane gives compactified 
phase portrait on the plane. In practice, polar or projective coordinates are introduced on the plane. Due to this 
central projection which was introduced by Poincare we can observe how critical points representing asymptotic states 
of the system are spread out along the equator. Moreover if r.h.s. of the system are given in polynomial form then 
we have simple test of structural stability of the system. In particular, a vector field on the Poincare sphere will 
be structurally unstable if there are non-hyperbolic critical points at infinity or if there is a trajectory connecting 
a saddles on the equator of Poincare sphere. In opposite case if number of critical points and limit cycles is finite 
dynamical system is structurally stable and therefore generic in the Peixoto sense. 

Introducing new variables (r, 9) with compactification at infinity: 

r 

ip — cos 6*, 

1 — r 
r 

y =-j sine, 

1 — r 
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1.5 



.3 




-0.3 



FIG. 1: The phase portrait for ^ — 3/50 < 3/25, it represents general behavior around the critical point of a stable node type, 
WaH = — 1 at this point. 

we receive dynamical system in the form 

61' = - i cos^ 6'|r2 sin^ 6'((1 - + 6^{1 ~ cos^ 9^ + 

+ sin 61 cos 61 (^(1 - rf + (1 - 6C)(6Cr^ cos6l(cos6' + sin6') + sin^ 6)^ + 

+ (^(1 - r)2 + (1 - sin2 + 6^r2(sin6l + cos 61)2^ (^2r2 sin6' cos 6* - (1 - rf +6^r^ cos6l(cos6l - sin6')) |, 
r' ^^r{l - r) sin 61 cos eijr^ cos^6l(^(l - r)^ + 6^(1 - 6^?-^ cos^ 6*)- 

-r^ sin 61 cos 61 (^(1 - + (1 - 6^){6^r^ cos 61 (cos 6* + sin 6') + sin=^ 6)^- 

- (^(1 - rf + (1 - 60^^ sin^ ^ + 6Cr^(sin6l + cos6i)2) (2r^ sin 6* cos - (1 - rf + 6^r^ cos6'(cos6l - sin 6*)) |, 

where a prime denotes now differentiation with respect to d/dr = (1 — rfd/da (Figs. [U [71 [8]). 

The number and character of critical points at infinity depends on the value of ^. In the general case we have 
critical points (rg = 1, 9q) : 9 = ±7r/2 and solution to the equation 

6.^(7 - 24^) sin^ 6i cos 6i + 18.^ sin 9 cos^ 61 + 3(1-4^) sin'^ 9 + 36^^ cos^ 9 = 0. 

This equation can be easily solved for arbitrary coupling constant ^. On Fig. [9] we present location and number of 
critical point at infinity as a function of coupling constant ^. 



CONCLUSIONS 



In this paper we considered the non-minimally coupled phantom scalar field as the FRW model with dark energy 
driving the acceleration of the current universe. We called it the extended quintessence and formulated the FRW 
dynamics with this form of dark energy in terms of the autonomous dynamical system. Hence we study all evolutional 
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0.4 0-6 0-8 1 1.2 ,1.4. 



FIG. 2: Typical trajectories on the invariant 2-dimensional submanifold for the quadratic potential function and ^ = 17/75. 
The phase portrait represents the generic behavior around a focus type critical point. For 3/25 < (, < 1/6 we have the global 
attractor of this type in the phase space, for ^ > 1/6 the critical point of this type is no longer a global attractor. 




FIG. 3: The phase portrait for the degenerated case ^ — 3/25. 
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FIG. 4: Phase portraits with compactification at infinity for minimaiiy coupied phantom scalar field ^ = for various potential 
functions, a) U{tp) oc tp'^ , critical point at ?/>' = and ip = oo plays the role of global attractor, wx = — 1 at this point; b) 
?7(^/>) oc ip'" and a = 1; c) Ui^p) oc exp {—X'lp) and A = 1, critical point at = and = — oo is the global attractor of the 
dynamical system and wx ~ —1 at this point; critical point at i^' = 0, i/' = co is a saddle type; d) Ui^p) oc exp( — At/)^) and 
A = 1, ^ = 0, the critical point at finite domain at t/j ' = = plays the role of global attractor and wx = —1 at this point, the 
points at infinity -tp' ^ 0, ip ^ are a saddle type. 



paths of the extended quintessence in the phase space for all admissible initial conditions. The structure of the phase 
space crucially depends on the constant of the non-minimal coupling because right-hand sides of the dynamical system 
depends on this parameter. We investigate bifurcations under changing of this parameter to distinguish some generic 
evolutional routes to the FRW model with the cosmological constant. We found that dynamical system admits the 
invariant submanifold {ip, tp') and trajectories in the long time evolution approach to this submanifold of the system. 
The behavior of trajectories on this phase plane give us information how the FRW model with the cosmological 
constant appears as a final asymptotic state - a limit set. Our main conclusion is that there are principally two 
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FIG. 5: Dependence of xq on the coupling constant ^ for two first terms in the Taylor series of the equation of state parameters 
(|27p and (|35|) : light gray (red el.version) - first term, dark gray (blue el. version) - second term and two linear parametrizations 
in the scale factor (left panel) and in the redshift (right panel), values of wo and wi from [sj. The horizontal dashed line denotes 
unphysical limit on the values of xq < It is interesting that astronomical data together with assumption yo — — §a;o 

indicate the value of ^ parameter near conformal coupling ^ = 




FIG. 6: The phase portrait with compactification at infinity for ^ — 3/50. The critical point at infinity is a saddle type. 



generic evolutional scenarios of approaching the hmit set in the future: 

1. node type scenario, 

2. focus type scenario (repellor-attractor scenario). 

In the first type of scenario trajectories approach the critical point monotonically without oscillations. In the 
second type of scenario the approach to the final state is through the characteristic infinite number of damping quasi- 
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FIG. 7: The phase portrait with compactification at infinity for ^ = 3/20. The critical point of a focus type plays the role of a 
global attractor in the phase space. 

oscillations. Additionally in this case the final limit set is always a focus type and it is achieved by evolution from an 
unstable node. 

It is shown that the type of evolutional scenario describing the approach toward the stable limit set in the future 
crucially depends on the trace and the determinant of the linearization matrix of the system on the invariant subman- 
ifold. They can be simply expressed in the terms of slow-rolling parameters, i.e., in terms of geometry of a potential 
function of the phantom scalar field. It was also demonstrated how our conclusion depends on the choice of a special 
form of a potential function of the scalar field. We found that in the generic case two mentioned before scenarios are 
typical. We obtained the exact form of coefficient of the equation of state w{z) near the present epoch directly from 
the dynamics of the model. This gave us the form of w{z) which can be used in estimation of this parameter. Such a 
methodology is more appropriate in our belief than taking the ad-hoc form of w{z) usually in a linear form. 
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